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ABSTRACT 


Filmwise  condensation  of  steam  on  externally-finned  tubes 
is  a  very  complex  process.   Recent  experiments  have  shown 
that  enhancement  ratios  (ratio  of  steam-side  heat-transfer 
coefficient  to  that  of  a  smooth  tube  having  the  same  diameter) 
exceeded  the  area  enhancement  produced  by  the  fins.   Moreover, 
the  enhancement  ratios  for  fully  flooded  tubes  exceeded  the 
values  predicted  by  a  simple,  one-dimensional  conduction 
model  by  a  factor  of  2  to  4 .   A  new  two-dimensional  conduction 
model  was  developed,  which  showed  that  the  one-dimensional 
model  overpredicted  the  two-dimensional  results  for  high 
conductivity  tube-metals  such  as  copper  by  as  much  as  13%. 
The  two-dimensional  model  also  showed  that  variations  in  fin 
thickness  or  spacing  can  result  in  an  overprediction  by  the 
one-dimensional  model  of  the  two-dimensional  results  by  as 
much  as  twenty- two  percent. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Because  the  size  and  weight  of  main  propulsion  and 
auxiliary  systems  are  inversely  proportional  to  the  ease  in 
which  heat  or  energy  is  exchanged,  a  major  thrust  of  the 
research  and  development  program  of  the  U.S.  Navy  is  toward 
the  design  and  development  of  equipment  that  perform  at 
higher  efficiency  with  a  reduction  in  size  and  weight.   In 
particular,  one  area  of  great  interest  is  the  reduction  in 
size  and  weight  of  steam  condensers. 

Current  equipment  is  of  the  bare  or  smooth  tube  type  of 
heat  exchanger  in  which  the  steam  is  condensed  on  the  outside 
of  the  tubes  (the  shell  side) .   The  heat-transfer  effective- 
ness is  governed  by  the  amount  of  surface  and  the  overall 
resistance  to  the  flow  of  heat.   This  overall  resistance,  in 
turn,  is  the  sum  of  the  vapor-side,  tube  wall,  and  water-side 
thermal  resistances.   Reduction  of  any  or  all  of  these 
resistances  results  in  an  increase  in  heat-transfer  performance 
with  an  associated  reduction  in  condenser  surface  area  which 
translates  into  decreased  size. 

Both  the  vapor-side  and  tube-side  thermal  resistances  can 
be  considered  as  the  reciprocal  of  the  so-called  hS  product 
where  h  is  the  coefficient  of  heat  transfer  and  S  is  the 
surface  area  through  which  the  heat  flows.   Because  of  the 
curved  surface  of  the  exchanger  tubing,  the  vapor-"space" 
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and  tube-side  surfaces  are  not  equal  and  an  adjustment  must 
be  made  to  the  tube-side  heat-transfer  coefficient  to  allow 
for  a  design  on  the  basis  of  the  vapor  or  shell-side  surface. 

It  is  a  fact  that  the  vapor-side  hS  product  (reciprocal 
of  resistance)  can  be  as  high  as  two  or  three  times  the  tube- 
side  hS  product.   Thus,  it  would  appear  reasonable,  in  an 
attempt  to  decrease  the  overall  resistance,  to  attack  and  try 
to  increase  the  hS  product  on  the  tube  side.   Unfortunately, 
little  can  be  done  in  the  arena  of  increasing  h  significantly 
on  either  side  although  the  vapor  side  has  greater  possibili- 
ties in  this  regard  and,  by  its  very  location,  not  much  can 
be  done  with  regard  to  the  increase  of  the  tube-side  surface. 

It  is  for  the  foregoing  reasons  that  recent  efforts  have 
focused  on  the  use  of  externally-finned  tubes  in  steam  con- 
densers to  increase  the  flow  of  heat  by  increasing  the  vapor- 
side  hS  product  [Refs.  1,2,3]. 

For  the  past  few  decades,  it  has  been  accepted  that  the 
use  of  externally-finned  tubes  would  be  inappropriate  mainly 
because  of  the  liquid  retention  (flooding)  that  occurs  between 
fins  (Figure  1.1).   This  flooding  is  a  result  of  the  relatively 
high  surface  tension  of  water  (at  50°C  water  has  a  surface 
tension  of  0.068  N/m,  while  refrigerants  possess  values  about 
0.015  N/m).   One  of  the  unavoidable  results  of  flooding  is  the 
considerably  lower  heat  transfer  through  the  flooded  portion 
than  that  through  the  unflooded  portion.   For  the  case  of  an 
externally  finned  copper  tube  with  steam  as  the  condensate, 

14 


External    Diameter   of    Fin^ 


Root    Dia 


Tube  W 


Retained  Condensate 


,     Retention  Angle  (l//) 


■igure  1.1.   Schematic  of  Condensate  Retention  on  Finned  Tubes 
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the  copper  has  a  thermal  conductivity  600  times  greater  than 
that  of  water  (38  5  W/m«K  for  copper  compared  with  0.6  W/m-K 
for  water) .   Due  to  its  overall  importance  in  the  heat- 
transfer  scheme,  the  flooding  of  the  lower  portion  of  the  tube 
by  condensate  has  been  studird  theoretically  and  experimentally 
by  several  investigators  [Refs.  4,5,6,7,8].  These  studies  have 
shown  that  the  flooding  angle  depends  mainly  on  the  fin 
spacing  and  the  ratio  of  surface  tension  to  density  of 
condensate  as  given  by  the  following  equation  for  the  case  of 
rectangular-section  fins  [Ref .  5] : 


H>     =     cos"1   [1 —-]  (1.1) 

r  pgDs 


Recent  work  by  Yau  et  al.  [Ref.  9]  and  Wanniarachchi  et  al.  [Ref.  1] 
has  shown  that  even  in  the  cases  with  considerable  flooding, 
the  heat-transfer  performance  exceeds  that  expected  on  the 
basis  of  area  increase  due  to  the  fins. 

In  1948,  Beatty  and  Katz  [Ref.  10]  developed  the  first  model 
to  predict  the  condensation  heat-transfer  coefficient  for 
horizontal,  externally-finned  tubes.   Their  model  ignores 
the  existence  of  flooding.   Although  their  model  agrees  to 
within  10  percent  for  condensing  fluids  possessing  low  surface 
tensions,  it  tends  to  be  rather  inaccurate  for  high-surface- 
tension  fluids,  such  as  water.   Several  other  investigators, 
Rudy  and  Webb  [Ref.  2] and  Owen  et  al.  [Ref.  8]  attempted  refinement 
of  the  Beatty  and  Katz  model  but  still  were  unable  to  predict 
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the  heat-transfer  coefficient  to  a  reasonable  degree  of 
accuracy. 

Rudy  and  Webb  used  the  Beatty  and  Katz  model  for  the 
unflooded  portion,  while  neglecting  heat  transfer  through  the 
flooded  portion.   With  this  modification,  Rudy  and  Webb  under- 
predicted  data  for  low-surface-tension  fluids  by  up  to 
50  percent.   Owen  et  al .  attempted  to  improve  the  Rudy  and 
Webb  model  by  allowing  heat  transfer  through  the  flooded 
portion.   In  order  to  compute  heat  transfer  through  the 
flooded  portion,  they  assumed  parallel  conduction  paths 
through  the  fins  and  the  condensate  that  is  flooded  in  the 
interfin  space  as  shown  in  Figure  1.2.   With  this  modification, 
they  improved  the  predictions  slightly;  underprediction  for 
low-surface-tension  fluids  decreased  to  about  4  0  percent. 
However,  as  shown  by  Honda  and  Nozu  [Ref.  11], both  of  these 
models  overpredicted  the  steam  data  of  Yau  et  al.  by  a  factor 
of  up  to  3.   The  fact  is  that  the  use  of  Beatty  and  Katz 
model,  which  does  not  take  surface  tension  forces  into  consi- 
deration for  the  unflooded  portion,  may  be  responsible  for 
this  large  discrepancy. 

The  most  recent  and  accurate  attempt  at  modelling  the 
problem  was  accomplished  by  Honda  and  Nozu  [Ref.  11]  .  Their 
resulting  solution,  a  fourth-order  differential  equation  for 
the  film  thickness,  was  solved  numerically  subject  to  reasonable 
boundary  conditions.   An  agreement  to  within  ±20%  was  achieved 
by  their  model  in  predicting  data  experimentally  obtained 
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Fiaure  1.2.   One-Dimensional  Parallel  Conduction  Kodel 
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for  11  fluids  and  22  tubes.   However,  the  disagreement  with 
steam  condensation  data  was  as  high  as  40%. 

To  understand  the  condensing  heat  transfer  performance 
more  clearly,  it  is  convenient  to  define  the  enhancement  ratio 
as  proposed  by  Wanniarachchi  et  al.  [Refs. 1,2].  The  enhancement 
ratio  is  defined  as  the  ratio  of  the  steam-side  heat-transfer 
coefficient  of  a  finned  tube  to  that  of  a  smooth  tube  with  the 
same  diameter  as  the  finned  tube  root  diameter  at  the  same 
heat  flux. 

The  present  investigator  has  been  mainly  motivated  by  the 
heat  transfer  through  a  finned  tube  that  is  completely  flooded 
by  condensate.   As  shown  by  Winniarachchi  et  al.  [Ref.  2],  a  finned 
tube  with  a  fin  spacing  of  0.5  mm,  and  a  height  and  thickness 
of  1.0  mm  and  1.0  mm  respectively  showed  a  considerable 
improvement  in  heat-transfer  performance  even  though  this  tube 
was  completely  flooded  by  condensate.   In  fact,  this  tube 
showed  an  enhancement  ratio  of  about  2.5  under  vacuum  (85  mmHg) 
and  4.5  at  atmospheric  pressure.   These  ratios  even  exceeded 
the  area  enhancement  (Af/A  ) ,  which  was  2.11.   If  one- 
dimensional  conduction  is  assumed  (see  Owen  [Ref .  8] ) ,  the  expected 
enhancement  ratio  was  computed  to  be  about  0.8.   This  unex- 
pectedly high  enhancement  ratio  may  be,  at  least  in  part, 
explained  by  two-dimensional  conduction  and  convection  in 
the  interfin  space. 
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B.   OBJECTIVES 

The  main  objectives  of  this  thesis  are: 

1.  Develop  a  two-dimensional  numerical  model  of  a  horizontal- 
finned  tube  fully  flooded  with  steam  condensate.   The 
program  using  the  model  shall  be  interactive  to  allow 
different  cases  of  fin  thickness,  fin  spacing,  and 
tube-metal  thermal  conductivity  to  be  studied, 

2.  Using  a  thermal  analyzer  computer  program,  determine 
the  temperature  distribution  plots  and  heat  flux  for 
each  case, 

3.  Compute  enhancement  ratios  for  each  one  of  the  cases 
studied  and  compare  these  ratios  with  those  computed 
by  using  a  one-dimensional  model, 

4.  Compare  the  one  and  two-dimensional  enhancement  ratios 
and  comment  on  the  importance  of  the  role  of  two- 
dimensional  heat  transfer  in  the  overall  heat-transfer 
performance, 

5.  Determine  the  value  and  usefulness  of  the  two-dimensional 
model  as  an  aid  in  the  prediction  of  enhancement  ratios 
for  various  cases. 
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II.   DEVELOPMENT  OF  THE  TWO  DIMENSIONAL  MODEL 

A.   THE  PRINCIPLE  OF  THE  FINITE  DIFFERENCE 

1.   Establishment  of  the  Spatial  Increment 

Consider  the  three-dimensional  Laplace  steady-state 
equation  for  the  spatial  temperature  distribution  in  a 
configuration 


2      2      2 
L_|  +  9T  +  3T   =   Q  (2>1) 

8x     8y     3z 


It  may  be  solved  analytically,  subject  to  certain  boundary 
conditions,  and  the  solution  is  called  an  analytical  solution 
The  finite-difference  solution  is  arrived  at  numerically  and 
applies  to  a  network  of  points  called  nodes,  with  the 
temperatures  at  the  nodes  being  estimated.   The  numerical 
solutions  obtained  by  analytical  and  finite-difference 
methods  can  be  made  to  agree  with  any  desired  accuracy. 

Consider  a  grid  of  nodes.   Applications  of  equation 
(2.1)  at  a  particular  node  whose  temperature  is  T   allows 
estimation  of  the  second  derivatives  in  terms  of  the 
temperature  T   and  the  six  related  temperatures  at  the 
adjacent  nodes  of  the  grid.   In  Figure  2.1  an  arbitrary 
function  T  =  f (x)  is  shown  which  considers  only  the  trace 
of  T  =  f 9 (x)  even  though  equation  (2.1)  shows  that  there 
will  also  be  a  functional  relationship  between  T  and  y, 
namely,  T  =  f  (y)  and  between  t  and  z,  namely,  T  =  f^(z) . 
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Figure  2.1.   Approximation  of  the  First  Derivative 
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If  y  and  z  are  held  constant,  the  first  derivative 
of  equation  (2.1)  may  be  approximated  by  the  change  in  T 
divided  by  the  change  in  x: 


8T   2    AT  +  AT       (T2  "  V  +  (To  '    V       T2  "  Tl 


8x      Ax..  +  Ax2  Ax,  +  Ax„  Ax,  +  Ax„ 

(2.2) 


An  estimate  of  the  second  derivative  is  obtained  by 
estimating  the  first  derivative  in  each  of  the  two  intervals 
and  dividing  the  difference  of  the  first  derivatives  by  the 
total  interval  over  which  the  first  derivative  estimates 
apply.   The  total  interval  will  be  the  distance  between  the 
centers  of  the  two  intervals,  Ax,  and  Ax„ .   Hence 


32T      AT2/Ax2  -  AT1/Ax1 

2   ~     I        1 
8x         j  Ax2  +  j   Ax, 


(2.3) 


If  Ax,  =  Ax2  =  Ax, 


92T      (T2-TQ)/Ax  -  (Tq-T1)/Ax      Ti  +  T2  "  2To 

"  2   ~  Ax  ™         Ax 

3x 


A  similar  development  using  T  =  f 2 (y)  with  temperatures 
T-,  and  T.  yields 


a2T      T 3  +  T   -  2T 

i_4   =   _^ 5 °  (2.4) 

2  Av 
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and  for  T  =  f->(z)  with  temperatures  TV  and  T, 


^T  T5  +  T6  -  2To 

3z2   ~         iz  ' 


so  that  equation  (2.1)  may  be  approximated  by 


T,  +  T~  -  2T     T0  +  T„  -  2T     Tc  +  Tr  -  2T 

* + 7 °  + 1 °  "  0  (2-6) 

Ax  Ay  Az 


and  if  Ax  =  Ay  =  Az,  the  result  is 


T,  +  TV  +  T,  +  TV  +  TV  +  TV  -  6T    =   0 
12     3     4     5     6      o 


Further  justification  of  the  result  indicated  by 
equation  (2.3)  comes  from  a  consideration  of  a  Taylor  Series 
Expansion.   Recall  that  the  Taylor  Series  Expansion  is  given 
by 


f(x)  =  fif(x-a)0  +  ^j^(x-a)  +  i^(x-a)2 

+  i^i(x-a)3  +  ...  +  ^l(x-a)n  (2.7) 


Two  cases  will  be  examined.   First  for  f (x)  =  T(x) 
at  x  =  x, ,  then 

T(x=x1)  =  T1 

T(x=x  )  =  T 
O      O 
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a  =  x 


x-a  =  Ax 


and  with  these  in  equation  (2.7) 


dT 
1    o    dx 
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a  2      1  d3T 

AX        +  -pr 
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6  A       3 

dx 


Ax" 
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Ax    +  .  .  . 


x 


o 


(2.3) 


Then  for  f(x)  =  T(x)  at  x  =  x. 


T(x=x2)  =  T2 

T(x=x  )  =  T 
o      o 


a  =  x 

x-a  =  -Ax 

With  these  in  equation  (2.7) 


T„  =  t   -  ~ 


d'r 
dx 


1  d2T 

Ax   +   tt 


x 


2  dx2 


Ax 


1  d3T 


x 


6,3 

dx 


Ax" 


x 


o 


4 
1  d  T 


24  ,  4 
dx 


Ax 


(2.9) 


x 


When  equations  (2.8)  and  (2.9)  are  added,  one  obtains 


2  4 

d  T     2     1  d  T     4 

T2  +  Tl  =  2To  +  TT(ix»  +  n   7T(Ax)   + 

dx  dx 


(2.10) 
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d2T 

ax2 

= 

Tl   +  T2   -   2To 

xo 

Ax2 

If  Ax  is  small,  then  all  terms  above  the  second  order 
in  Eq.  (2.9)  may  be  considered  negligible.   Thus  upon 
rearrangement,  it  is  seen  that 


(2.11) 


which  is  a  result  that  is  identical  to  eauation  (2.3).   Of 
course  a  similar  development  will  provide  a  substantiation 
of  equations  (2.4),  (2.5)  and  (2.6). 
2 .   Node  to  Datum  .Analysis 

Using  the  approximations  for  the  second  derivative, 
a  node  analysis  can  be  performed  to  predict  the  temperature 
distribution  in  a  mathematical  model  of  a  physical  configura- 
tion.  With  a  node  analysis,  the  configuration  of  interest 
is  divided  into  small  but  finite  subvolumes  with  each  sub- 
volume  being  considered  as  isothermal.   It  may  be  then 
assumed  that  the  geometric  center  of  each  subvolume  is 
representative  of  that  subvolume  and  this  geometric  center 
is  called  the  node-point  or  node. 

A  node  analysis  consists  of  an  application  of  the 
First  Law  of  Thermodynamics  to  generate  energy  balances  for 
each  node.   This  is  done  for  each  node  by  equating 

HEAT  IN  +  HEAT  GENERATED  =  HEAT  OUT  +  HEAT  ACCUMULATED    (2.12) 


26 


For  the  two  dimensional  case  displayed  in  Figure  2.2,  the 
development  of  the  node  equations  for  node  point  (j,k)  is 
shown  in  Figure  2.3.   Assuming  heat  flow  by  conduction  and 
a  nodal  dissipation,  q,  a  steady  state  or  First  Law  energy 
balance  requires  that  the  sum  of  all  heat  flows  (Watts  or 
BTU/hr)  leaving  the  node  equals  zero.   Thus  with  a  heat  input 
as  a  negative  heat  flow  outward, 

qx  +  q2  +  q3  +  q4  -  ^i  =  °         (2.13) 

In  terms  of  the  temperatures  in  Figure  2.4  and  in  terms  of  the 
Fourier  Law 


q   =   M^T).  (2,14) 


«i    =    ^*<VTi> 


<J2  "  T5F"V*2> 


a   =  MX  5(T  _T  ) 
q3      Ax  °  uo  i3) 


^4  =  T7  «(VT4> 


and  when  these  are  inserted  into  equation  (2.13),  the  result 
is 
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Figure    2.2.       Two-Dimensional   Case 
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Figure  2.3.   Node  Equation  Point  (j,k 
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Figure    2.4.      Node   Equation 
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^*<VV  +T7*<VT2>  +  ^<W 


+   ^  5<VV       =      <3i  <2-15' 


If  two  conductances  are  defined  by  Kl  =  KAy6/Ax  and 
K2  =  KAX6/AY  then 

KICTq-T^)  +  K2(TQ-T2)  +  K1(T  -T3)  +  K2(Tq-T4)  =  q. 

(2.16) 

and  if  there  is  no  heat  generation  (q.=0),  the  equation 
becomes 


Tq(2K1  +  2K2)  -  KK^)  -  K2(T2)  -  Kl(T3)  -  K2(T4)  =  0 

(2.17) 


For  the  case  where  there  is  heat  flow  by  both 
conduction  and  convection  as  in  Figure  2.5,  with  a  heat 
input  as  a  negative  flow  outward, 

ql  +  q2  +  q3  +  qc  "  qi   =   °  (2.18) 

For  the  conduction  terms,  q, ,  q^  and  q_,  equation  (2.14)  is 
still  applicable.   However,  for  the  convection  at  the  boundary 
as  indicated  by  q   (Figure  2.5) 


qc   =   hAy5(AT)  (2.19) 
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Figure  2.5.   Node  on  Convecting  Boundary 
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In  terms  of  the  temperatures  in  Figure  2.5 


K(Ax/2)  f  ,m   m 
ql   =   W^   6  (To"Tl 


q9   =   ^  6(T  -T- ) 
^2       Ax     o   2 


q    =   K(Ax/2      _ 

Ay   °  uo  x3; 


qc   =   hAy  6(To-Ts) 


and  when  these  are  inserted  into  equation  (2.18),  the  result 
is 


K(Ay2>  6  (T  -T.  )  +  m   «  (Tn-T, )  +  K(A*/2)  5  (T  -T,) 
Ay      o   1      Ax     0   2        Ay       o   3 


+  hAy6 (T  -T  )   =   q.  (2.20) 

-'os        1 


Letting  Kl  =  K(^/2)  ,  K2  =  ^|^  and  K3  =  hAy6  ,  then 


K1(.T  -T,)  +  K2(TQ-T2)  +  Kl(TQ-T3)  +  K3(Tq-Ts)   =   qi 

(2.21) 
which  can  be  simplified  to 


(2K1  +  K2  +  K3)T   -  KIT,  -  K2T2  -  KlT3   =   qi      (2.22) 
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and  once  again  if  there  is  no  heat  generation  (q.  =0),  the 
equation  becomes 

(2K1  +  K2  +  K3)T   -  KIT,  -  K3T„  -  KIT   -  K3T    =   0      (2.23) 

For  the  case  of  a  node  on  an  adiabatic  boundary  with  no 
heat  generation  as  in  Figure  2.6 

qx  +  q2  +  q3  =  °  (2.24) 

In  terms  of  the  temperature  in  Figure  2.6  and  in  terms  of  the 
Fourier  Law 


qi   =  KUx^^, 
^2   =   T^'VV 


=   K(Ax/2) 
^3         Ay       o   3 


and  when  these  are  inserted  into  equation  (2.24),  the  result 
is 


K(Ax/2) 


*<VT1>  +  ^«VV  +H^^'VT3» 


Ay      o   1'     Ax     o   2'       Ay 

(2.25) 

r  4-4--    vi         K(Ax/2)  .    ,  „_    KAy6 
Letting  Kl  =  — 5 — jr — -  6  and  K2  =  —r1 — 

^  Ay  Ax 


(2K1  +  K2)TQ  -KlT^  -  K2T2  -KIT    =   0  (2.26) 


It  is  easy  to  see  that  when  one  proceeds  in  this 
manner  for  all  nodes,  an  nxn  set  of  simultaneous  algebraic 
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Figure  2.6.   Node  on  Adiabatic  Boundary 
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equations  will  result.   These  may  be  represented  in  matrix 
form  as 

[K] [T]   =   [Q]  (2.271 

where  [K]  is  an  nxn  matrix  of  coefficients  and  IT]  and  [Q] 
are  nxl  column  vectors  representing  the  node  temperatures 
and  heat  inputs  respectively. 

The  solution  of  equation  (2.27)  is  handled  by  the 
computer  program  TVSSI  [Ref.  12:Appendix  A] . 

B.   MODEL  BUILDER  TVIN 

In  order  to  provide  input  to  program  TVSSI,  it  was 
necessary  to  create  a  model  of  the  configuration  to  be 
analyzed.   The  overall  model  of  the  condensing  region 
incorporated  three  major  features.   Because  the  system  to 
be  modelled  was  an  externally  finned  tube  with  the  external 
space  flooded  by  condensate,  the  azimuthal  portion  of  the 
model  was  chosen  to  be  1/2tt  radians.   This  is  shown  in 
Figure  2.7  and  because  the  areas  of  each  node  contain  a  tt 
term,  the  use  of  1/2tt  simplified  the  construction  of  the  node 
equations  by  permitting  cancellation  of  the  tt  terms.   The 
second  feature  required  that  the  maximum  number  of  nodes  to 
be  used  for  the  model  was  300,  due  to  input  restrictions 
(the  dimension  statement)  of  the  main  thermal  analyzer.   The 
final  requirement  necessitated  that  the  program  be  written  so 
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Figure    2.7.      Model   Azimuthal    Section 
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as  to  allow  the  user  to  input  the  variables  listed  in  Table  I 
interactively . 

The  basic  model  is  shown  in  Figure  2.8.   The  upper 
horizontal  edge  represents  the  condensate/vapor  interface 
while  the  lower  horizontal  edge  represents  the  inner  wall  of 
the  tube.   Heat-transfer  coeffients  were  specified  for  the 
inner-tube  wall  and  the  vapor-condensate  interface.   The 
model  assumes  that  the  heat  is  transferred  across  the  retained 
condensate  by  conduction  only.   Due  to  symmetry,  the  dotted- 
dashed  lines  at  the  two  vertical  edges  of  the  figure  represent 
an  isolation  of  a  fin-condensate-tube  entity  into  a  repeating 
section.   These  lines,  of  course,  represent  adiabatic 
boundaries.   Thus,  the  model  assumes  that  the  only  heat 
transferred  into  or  out  of  the  model  occurs  at  either  the 
condensate/vapor  interface  or  the  inner  tube  wall. 

The  first  model  constructed  consisted  of  299  nodes  with 
13  columns  and  23  rows.   The  outermost  rows  and  columns 
consisted  of  half  node  subvolumes.   This  facilitated  the 
placement  of  the  node  centers  at  both  the  inner  tube  wall 
and  upper  condensate/vapor  interface  in  addition  to  the 
boundary  of  the  outer  tube  wall  and  fin  perimeter.   Nine  rows 
were  placed  in  the  condensate  region  above  the  fin  and  in  the 
region  including  the  fin  and  adjacent  condensate.   The 
remaining  four  rows  were  placed  in  the  tube  wall  region. 
In  addition  the  fin  region  contained  4  columns  of  nodes 
while  the  remaining  eight  columns  were  placed  in  the  inter-fin 
space. 
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Figure    2.8.       Basic   Model 
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TABLE  I 
INPUT  VARIABLES  FOR  MODEL  BUILDER  WIN 


T.     -  Temperature  of  fluid  inside  tube 


T     -  Temperature  outside  condensate  layer 


r.     -  Inner  radius  of  tube 

1 


r     -  Outer  radius  of  tube 
o 


h.     -  Heat  transfer  coef .  at  inside  tube  wall 

1 


h     -  Heat  transfer  coef.  at  condensate  layer 


k      -  Thermal  conductivity  at  tube  and  fin 


k,  -  Thermal  conductivity  of  condensing  liquid 

e  -  Fin  height 

t  -  Fin  width 

s  -  Spacing  between  fins 

z  -  Height  of  condensate  above  tip  of  fin 
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The  location  of  these  columns  and  rows  determined  the 
height  and  width  increments  of  the  various  subvolumes, 
depending  on  their  location  in  the  model.   For  example,  if 
the  tube  wall  thickness  is  four  millimeters,  and  the  tube 
wall  area  contains  four  rows,  the  height  increment  for  a 
node  in  this  region  would  be  4  mm/4  rows  =  1  mm.   The  width 
increments  for  a  node  would  be  computed  in  identical  fashion 
by  dividing  the  size  of  the  particular  region  by  the  number 
of  columns  occupying  that  region.   Figure  2.9  identifies  all 
of  the  different  height  and  width  increments  utilized  by  the 
model. 

C.   EVALUATION  OF  RESULTS 

1 .   Verification  of  Accuracy 

Verification  of  the  model  consisted  of  using  the 
temperatures  generated  by  TVSSI  (Using  TVIN  as  input)  and 
the  conductance  values  employed  in  TVIN  to  compute  the  heat 
flow  across  three  boundaries  in  the  model : 

(1)  condensate  upper  surface 

(2)  fin  perimeter  and  tube  outer  wall 

(3)  tube  inner  wall 

The  foregoing  heat  flows  were  computed  in  each  case 
by  using  the  node  temperatures  generated  by  TVSSI  with  the 
associated  internode  conductance  values  used  in  the  input 
generated  by  TVIN.   These  heat  flows,  of  course,  must  be 
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Figure  2.9.   Model  Increments.   (A)  Height  and  (B)  Width 
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equal  and  it  was  observed  at  the  outset  that  13x23  model 
yielded  certain  unacceptable  discrepancies  in  the  three  values 

Because  the  13x23,  299  node  model  was  found  to  be 
inaccurate,  it  was  determined  that,  in  order  to  achieve  the 
requisite  accuracy,  slight  adjustment  of  the  number  of  rows 
and  columns  was  required  from  run-to-run. 

TVIN  was  reprogramed  to  allow  flexibility  in  the 
number  of  rows  and  columns  utilized  in  the  model  as  well  as 
placement  of  the  nodes  in  the  fin,  tube  or  liquid  regions. 
The  first  change  instituted  was  to  fix  the  number  of  rows 
contained  in  the  condensate  layer  above  the  fin  tip  at  two 
rows.   This  action  was  acceptable  because  the  original  height 
of  this  layer  was  only  0.0  8  mm  compared  with  a  1.0  mm  fin 
height  and  3.175  mm  tube  wall  thickness.   Fixing  the  number 
of  rows  at  two  permitted  the  use  of  additional  rows  at  the 
fin  and  tube  levels.   Program  TVIN  was  changed  to  inter- 
actively accept  a  specified  number  of  columns  (see  Appendix 
A) .   The  number  of  rows  was  obtained  by  making  the  product 
of  rows  and  columns  as  close  as  possible  to  300. 

A  description  of  the  menu-driven  procedure  for  the 
selection  and  allocation  of  rows  and  columns  for  a  296-node 
model  example  now  follows: 

If  the  response  to  the  query  'How  many  columns  are 
desired?'  is  8,  the  program  will  provide  300/8  =  37.5  or  37 
rows.   This  yields  an  8x37  =  296  node  model. 

Suppose  the  response  to  the  query  'How  many  rows  are 
desired  for  the  tube  wall?'  is  5.   Then  the  program  will 
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allocate  37-5  =  32  rows  for  the  fin  and  condensate  layer 
above  the  fin. 

However,  the  number  of  rows  in  the  condensate  layer 
above  the  fin  tip  is  fixed  at  2.   Thus  the  number  of  rows  in 
the  fin  becomes  32-2  =  30.   The  result  is  an  8x37,  296  node  model 
A  display  of  the  foregoing  model  set-up  is  shown  in 
Figure  2.10. 

2.   Model  Verifier  Ql 

As  mentioned  in  the  previous  section,  once  an  input 
file  had  been  created  by  TVIN,  and  the  node  temperatures 
were  calculated  by  the  thermal  analyzer,  TVSSI,  the  solution 
was  verified  for  correctness  and  accuracy.   This  check 
consisted  merely  of  a  comparison  of  the  heat  flow  (watts), 
across  three  surfaces  in  the  model.   An  error  of  3.5%  was 
established  as  the  maximum  allowable  for  a  solution  to  be 
considered  acceptable.   The  model  verifier  Program  Ql  (see 
Appendix  B) ,  is  an  offshoot  of  the  model  builder  TVIN  as  it 
interactively  accepts  the  same  input  data  required  to  build 
the  model.   Once  all  input  parameters  are  set,  the  program 
generates  and  stores  each  inter-connecting  node  conductance 
value. 

In  the  next  step,  all  node  temperatures  computed  by 
the  thermal  analyzer  are  read  and  stored,  by  location  in  the 
model.   The  node  numbers  of  the  condensate  outer  layer,  fin 
perimeter  and  outer  tube  wall  and  the  inner  tube  surface  are 
determined  for  each  model  via  an  algorithm  which  requires  as 
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Figure  2.10.   8x37,  296  Node  Model 
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input,  the  number  of  columns  and  rows  being  used  in  the  model 
along  with  their  placement.   In  addition  to  determining  the 
nodes  on  these  surfaces,  the  algorithm  also  computes  the 
node  number  of  each  of  the  adjoining  nodes  in  the  direction 
of  the  heat  flow.   Computation  of  the  heat  flow  by  program  Ql 
is  then  simply  a  matter  of  taking  each  surface  node  conductance 
and  multiplying  it  by  the  temperature  difference  between  the 
surface  node  and  its  adjoining  node.   The  summation  of  all 
the  surface  node  values  results  in  the  total  heat  flow  across 
the  particular  surface. 

D.   THERMAL  ANALYZER  TVSSI 

The  computer  program  TVSSI  used  in  obtaining  the  temperature 
distributions  was  adapted  from  program  TVSS2  listed  by  Kern 
and  Kraus  [Ref.  13: Appendix  2a].   This  adaptation  consisted 
of  modifying  TVSS2  to  compute  solutions  in  the  SI  system  and 
accept  an  input  file  specifically  designed  by  model  builder 
TVIN  for  each  externally  finned  tube  case.   This  program  also 
had  to  be  adapted  to  the  IBM  30  33  AP  system. 

The  program  utilizes  the  Cholesky  decomposition  scheme, 
and  because  of  the  linearization  of  the  radiation  terms  (not 
used  in  this  study),  the  program  is  iterative.   A  more  indepth 
description  of  the  Cholesky  Scheme  pertaining  to  the  program 
is  discussed  by  O'Hare  [Ref.  12:pp.  48-52]. 

For  all  the  cases  which  were  studied,  the  following 
parameters  were  set  for  the  analyzer: 
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1.  The  initial  temperatures  at  which  all  nodes  were  set 
and  at  which  the  computer  began  the  iterative  process 
was  100°C. 

2.  An  eventual  accuracy  of  .05  between  the  final  and  next 
to  last  iterations  was  used. 

3.  A  radiation  coefficient  convergence  factor  Of  0.66667 
between  iterations  was  used.  This  is  required  by  the 
program  even  if  radiation  is  not  considered  as  a  mode 
of  heat  transfer  in  the  model. 

4.  The  maximum  number  of  iterations  that  the  computer  was 
allowed  to  perform  was  set  at  12  to  avoid  excessive 
computer  running  time. 

5.  A  damping  factor  of  0.8  was  set  as  an  initial  damping 
factor  based  on  the  number  of  non-linear  terms  in  all 
of  the  node  equation. 

The  ability  of  TVSSI  to  correctly  solve  a  numerical  model 

was  verified  by  providing  as  input  a  data  set  listed  by  Kern 

and  Kraus  [Ref.  13:pp.  326-336]  for  an  extended  surface 

problem  where  the  solutions  were  derived  analytically.   The 

computer  generated  solution  agreed  completely  with  the 

analytic  solution  thereby  assuring  confidence  of  analyzer 

generated  temperatures. 
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III.   RESULTS  AND  DISCUSSION 

A.  INTRODUCTION 

During  this  thesis  effort,  15  different  cases  were  studied 
using  the  programs  TVIN  and  TVSSI.   Table  II  lists  the  param- 
eters that  remained  constant  in  all  cases.   A  summary  of  all 
cases  examined  and  their  results  are  listed  in  Table  III. 

B.  SMOOTH-TUBE  EVALUATION 

As  a  method  for  comparison  of  the  variation  in  heat- 
transfer  performance,  the  heat  flux  for  a  smooth  tube  with 
the  same  diameter  as  the  finned-tube  root  diameter  was  estab- 
lished as  a  baseline  value.   The  ratio  of  the  one-  and  two- 
dimensional  heat  fluxes  to  the  heat  flux  of  the  finned  tube 
with  fin  width  approaching  zero  were  designated  <j>.  and  $~ 
respectively.   The  equation  used  to  determine  the  finned-tube 
of  zero  fin  width  heat  flow  is  qiven  bv: 


(T  -T.) 
n  _  o   1 

y  D        D  (ln(D  /DJ)         -,  (3.1) 

o       o     o'  1    ,  e 


(hl>(V  +      2km       +  \+   h0 


C.   ONE  DIMENSIONAL  HEAT  FLOW 

The  one-dimensional  model  utilized  in  this  study  assumed, 
as  did  Owen  et  al.  [Ref.  8],  parallel  paths  for  heat  flow  across 
the  tube  wall,  fin  and  retained  condensate.   Furthermore,  a 
simplifying  assumption  was  made  that  heat  flow  is  purely  radial 
and  that  the  heat  is  transferred  across  the  retained  condensate 

4  8 


TABLE  II 
CONSTANT  PARAMETERS  FOR  MODEL  BUILDING  PROGRAM  TVIN 


h.   =   30000  W/m-K 


h    =   100000  W/m. K 
o 


r .   =   6.35  mm 

l 


r    =   9.525  mm 
o 


e    =   1.0  mm 


K,   =   0.6  W/m- K 


T    =   48.0  °C 
o 


T,   =   22.0  °C 
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by  conduction  only.   This  was  accomplished  by  modifying  pro- 
gram TVIN  in  such  a  way  so  that  the  data  file  it  created  con- 
tained only  those  internode  conductance  values  representing 
heat  flow  in  the  radial  direction  while  deleting  those  conduc- 
tance terms  pertaining  to  axial  heat  flow.   Programs  TVSSI  and 
Ql  were  utilized  to  compute  the  node  temperatures  and  heat 
flow  for  each  case.   Once  again,  a  maximum  error  of  3.5%  was 
used  when  heat  flows  generated  by  program  Ql  were  verified. 

D.  TWO  DIMENSIONAL  HEAT  FLOW 

The  heat  flow  for  the  two-dimensional  case  was  evaluated  by 
the  methods  discussed  in  Chapters  I  and  II. 

E.  EFFECTS  OF  FIN  SPACING  AND  FIN  WIDTH  ON  PERFORMANCE 

Heat  fluxes  were  computed  for  five  different  cases  of  fin 
width  and  spacing.   Three  cases  had  a  fin  width  of  1  mm  with 
spacings  of  0.5,  1.0  and  1.5  mm  respectively.   The  other  two 
cases  held  the  spacing  constant  at  0.5  mm,  while  the  fin  width 
was  adjusted  to  0.5  and  0.25  mm. 

1 .   Effects  of  Fin  Spacing 

This  section  presents  results  showing  the  variations 
of  heat  flux  having  fin  spacing  as  a  parameter.   Figures  3.1, 
3.2  and  3.3  present  results  for  tubes  having  a  fin  width  of  1.0 
mm  with  fin  spacings  of  0.5,  1.0  and  1.5  mm.   Each  figure  repre- 
sents the  same  spacing  variation  with  the  tube-metal  thermal 
conductivity  being  varied  from  385.0  W/m.K  for  copper,  through 
45.0  W/m  K  for  copper-nickel,  down  to  16.0  W/m-K  for  the  case  of 
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stainless  steel.   Similarly,  Figures  3.4,  3.5  and  3.6  present 
the  ratio  of  the  two-dimensional  heat  flux  to  the  one-dimensional 
heat  flux  (Q21) .   For  each  case  of  tube-metal  thermal  conduc- 
tivity, the  fin  spacing  of  0.5  mm  shows  the  best  heat-transfer 
performance,  while  the  worst  performance  occurs  at  a  spacing  of 
1.5  mm.   Figures  3.7,  3.8  and  3.9  present  isotherm  plots  for 
all  three  cases  of  thermal  conductivity  for  the  finned  tube 
have  a  0 . 5  mm  fin  spacing  and  1.0  mm  fin  width.   Moreover,  the 
ratio  Q21  has  its  smallest  value  at  a  find  spacing  of  0.5  mm 
and  increases  with  increasing  fin  spacing.   Enhancement  ratios 
for  the  various  cases  had  their  maximum  value  at  the  0 . 5  mm 
fin  spacing  and  this  ratio  showed  a  marked  decreaase  as  the  fin 
spacing  was  increased. 

The  results  indicate  several  important  facts.   The  first 
and  most  important  is  the  significant  role  the  two-dimensional 
heat  flow  plays  in  this  type  of  extended  surface  analysis.   The 
increasing  Q21  ratio  with  increasing  fin  spacing  shows  that  for 
the  cases  of  copper  and  copper-nickel  tube-metal,  the  axial 
heat  flow  into  the  fin  upper  side  surface  tends  to  raise  the 
temperature  of  the  fin  tip  resulting  in  a  choking  effect  of  the 
heat  flow  into  the  fin  tip  surface.   This  is  of  enormous  conse- 
quence because  in  each  and  every  case,  the  majority  of  the 
heat  absorbed  by  the  finned  tube  is  accomplished  at  the  fin 
tip.   Because  the  one-dimensional  case  allows  no  axial  heat 
flow  into  the  fin  side  surface,  a  greater  temperature  gradient 
occurs  at  the  fin  tip  resulting  in  a  greater  heat  flow.   As  the 
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Figure  3.7. 


Isotherm  Plot  For  Finned  Copper  Tube 
(Fin  and  Condensate  Regions) 
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Figure  3.8.   Isotherm  Plot  For  Finned  Copper-Nickel  Tube 
(Fin  and  Condensate  Regions) 
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Fiaure  3.9 


Isotherm  Plot  For  Finned  Stainless  Steel  Tube 
(Fin  and  Condensate  Regions) 
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fin  spacing  is  increased,  the  effect  of  the  fin  side  surface 
is  less  pronounced  since  the  poor  conductivity  of  the  condensate 
inhibits  axial  heat  flow  in  the  region  near  the  fin  space 
centerline.   This  results  in  a  more  closer  agreement  between 
the  one-dimensional  model  and  the  two-dimensional  version. 

This  conclusion  can  be  further  justified  by  comparing 
the  Q21  ratios  for  each  tube-metal  conductivity.   The  two- 
dimensional  effect  at  the  fin  tip  is  decreased  as  the  thermal 
conductivity  of  the  fin  material  decreases.   This  does  not 
mean,  however,  that  there  is  no  two-dimensional  effect  for  the 
cases  with  low  tube-metal  conductivity.   As  in  the  case  of 
copper,  the  two-dimensional  effect  tends  to  inhibit  the  heat 
flow  into  the  fin.   For  the  cases  of  low  metal  thermal  conduc- 
tivity, there  is  a  less  pronounced  two-dimensional  effect  at 
the  fin  tip,  but  there  is  a  significant  two-dimensional  heat 
flow  into  the  lower  portion  of  the  fin  side  surface  near  the 
tube  wall.   For  this  reason,  the  one-dimensional  model  tends 
to  underpredict  the  two-dimensional  model. 

As  seen  in  all  the  cases  for  the  copper-nickel  and 
stainless-steel  tube-metal  (except  the  copper-nickel  case 
with  0 . 5  mm  fin  spacing),  all  Q21  ratios  are  greater  than 
unity.   For  the  copper-nickel  case  of  0.5  mm  spacing,  the  reduced 
spacing  even  with  a  lower  metal  conductivity  of  45  W/m.K,  tends 
to  produce  a  two-dimensional  effect  at  the  upper  fin  tip. 
2.   Effects  of  Fin  Width 

The  results  of  the  variation  of  heat  flux  with  fin 
width  as  a  parameter  indicate  that  fin  width  plays  an  important 
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role  in  the  overall  heat-transfer  mechanism.   The  two  cases 
for  each  tube-metal  thermal  conductivity  investigated  were  fin 
widths  of  0.5  and  0.25  mm  respectively  with  fin  spacing  held 
constant  at  0 . 5  mm.   For  each  case  of  tube-metal  thermal  con- 
ductivity, the  fin  width  of  0.5  mm  displayed  one  of  the  highest 
heat  fluxes,  being  second  only  to  the  case  of  the  1.0  mm  fin 
width  with  0 . 5  mm  fin  spacing.   The  only  significant  two- 
dimensional  fin  tip  effect  can  only  be  observed  in  the  case 
of  the  copper  tube-metal  with  a  fin  width  of  0.5  mm. 

The  same  reasoning  as  in  the  previous  section  can  be 
applied  in  this  case.   Because  the  fin  tip  width  has  been 
reduced,  the  total  amount  of  heat  transferred  into  the  tip 
region  is  also  reduced  and  the  two-dimensional  fin  tip  choking 
effect  has  less  influence.   This  concurs  with  the  computed  Q21 
ratio  of  1.04.   By  reducing  the  fin  width  further  in  each  case 
to  0.25  mm,  each  of  the  resultant  heat  fluxes  is  less  than  the 
cases  with  0.5  mm  fin  width.   For  each  case  of  the  0.25  mm  fin 
width  the  Q21  ratio  approached  1.3,  signifying  the  comparatively 
large  two-dimensional  heat  flow  into  the  lower  fin  side  surface. 


63 


IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

1.  The  two-dimensional  numerical  model  of  a  horizontal 
finned  tube  fully  flooded  with  steam  condensate  developed  is 
a  useful  tool  in  the  study  of  the  relationship  between  one- 
and  two-dimensional  heat-transfer  effects. 

2.  Two-dimensional  heat  transfer  effects  play  a  major 
role  in  the  overall  heat  transfer  mechanism  in  the  case  of  an 
externally  finned  tube  fully  flooded  with  steam  condensate. 
It  has  been  determined  that  for  cases  of  high  tube-metal  con- 
ductivity, such  as  the  case  of  copper,  the  one-dimensional 
model  can  overpredict  the  results  of  the  two-dimensional 
model  by  as  much  as  13%,  while  variation  in  the  fin  width  or 
spacing  can  result  in  the  one-dimensional  model  underpredicting 
the  two-dimensional  result  by  as  much  as  twenty- two  percent. 

3.  Fin  width  and  fin  spacing  control  the  magnitude  of  the 
two-dimensional  effect  present  and  the  resulting  heat-transfer 
performance.   An  increasing  fin  width  to  fin  spacing  ratio 
will  correspond  to  an  increase  in  the  heat-transfer  performance. 

4.  The  one-dimensional  model  in  cases  of  high  tube-metal 
thermal  conductivity  exhibiting  a  fin  width  to  spacing  ratio 

of  greater  than  unity  will  overpredict  the  heat  flux  predictions 
of  the  two-dimensional  model. 

5.  Decreasing  the  tube-metal  thermal  conductivity  results 
in  a  degraded  heat-transfer  performance. 
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6.   The  two-dimensional  fin-tip  choking  effect  is  minimal 
in  low  thermal  conductivity  tube-metals  while  the  two- 
dimensional  heat  transfer  into  the  fin  lower  side  surface  is 
far  more  pronounced. 

B.   RECOMMENDATIONS 

1.  Use  the  two-dimensional  model  computed  enhancement  ratios 
for  the  various  cases  to  estimate  a  steam-side  heat-transfer 
coefficient  and  compare  with  those  methods  currently  being  used. 

2.  Using  the  model  builder  TVIN  and  thermal  analyzer 
TVSSI,  conduct  a  more  in-depth  study  of  the  effects  of  variation 
of  fin  width,  height  and  spacing  on  the  overall  heat-transfer 
performance. 

3.  Revise  the  two-dimensional  model  to  include  convection 
effects  in  the  retained  condensate  region. 
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AEP.EUDJLX.   A 
TVIN   COMPUTER    GENERATED     INPUT    ANALYZER    PROGRAM 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccc:cc 

C  PROGRAM    TVIN  C 

C  PROGRAM    T  VI M     IS   DESIGNED    TO    ACCEPT    VARIABLE  C 

C  DIMENSIONS    AND    INPJT    PARAMTERS    TO    CREATE  C 

C  AN    OUTPUT    DATA    FILE    OF    INTER_NODE    CONNECTIONS  C 

C  AND    CONDUCTANCE    TERMS.       BY    ADJUSTING    THE  C 

C  PRINT    STATEMENT,    THE    OUTPUT    WILL    EITHER  C 

C  REPRESENT    A    ONE    OR    TWO-DIMENSIONAL    MODEL.  C 

c  c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

c 

DIMENSION    CK!300,4),T(300) 

INTEGER    P(300,4) , CO ,N , X, ZW, ZW1, Z W2 ,Z W3 , VI , V2 ,W,U, 
DOUBLE    PRECISION    CTHM ,CT HL, HI ,HS ,RO, RI , A ,B, C , D, 
*VAL4,VAL14,VAL6,VAL16,VAL8, VAL18 , VAL5 , VAL7, V  ALIO , 
*A2,A3,A4,A5,WINC1 , WINC2 , H INC1 ,HI NC2, HI NC3 , AL INE , 
*TIN, TS,T,X1,X2,VAL2,VAL12,VAL20, Al ,RADH, RADL , EL 
CHARACTER    *79     TITLE 
CHARACTER    *79    FNAME 
C 

DATA     IOT,IN, IPR,I WR/6,5,4,8/ 
C 

OPEN(IPR,FILE=«PRN»  ,STATUS=  •NEW*  , FORM= • FORM ATTED'  , 
*IOSTAT=ICK) 
IFdCK    .NE.    OJ    WRITE(  IOT  ,20  ) 
C 

20         FORMAT! 'TROUBLE    OPENING    PRINTER    OUTPUT    FILE') 
C 

WRITEdOT,  17) 

17  FORMAT!/'     ENTER    THE    TITLE    OF    THIS    STUDY,    79     •, 

*  •    COLUMNS    MAX.       THIS    TITLE    WILL    APPEAR1,/,1     AT', 

*  •    THE    TOP    OF    EVERY    PRINTED    PAGE    OF    OUTPUT:1) 
READ(IN,18)    TITLE 

18  F0RMATU79) 
C 

WRITEdOT,  1) 

1  FORMAT!/'     ENTER    INNER    RADIUS    OF    TUBE:') 
READ(IN,2)    RI 

2  F0RMAT(BN,F10.5) 
C 

WRITEdOT, 3) 

3  FCRMAT(/«     ENTER    OUTER    RADIUS    OF    TUBE:') 
READ (IN, 2)    RD 

C 

WRITEdOT, 4) 

4  FORMAT!/'     ENTER    DISTANCE    FROM    OUTER    EOGE    OF     ', 

*  'FIN    TO    LIQUID    SURFACE:') 


READ(IN,2)     A 

WRITEdOT, 5) 

FORMAT!/'     EN1 
READ! IN, 2)    8 

WRITEdOT, 6) 

FORMAT!/'     EN" 
READ (IN, 2)    C 


*RITE!IQT,5) 

FORMAT!/'     ENTER    HEIGHT    OF    FIN') 


WRITEdOT, 6) 

FORMAT!/'     ENTER    DISTANCE    BETWEEN    FINS') 
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WRITE(I0Tt7) 

7  FORMATC/'     ENTER    WIDTH    CF     FIN') 
READCIN.2)    D 

WRITECIOT,8) 

8  FORMAT!/'     ENTER    HEAT    TRANSFER    COEFFICIENT    OF     \ 
*■     TUBE    INNER    SURFACE:  '  ) 

READUN,2)    HI 

WRITE(I0T,9) 

9  FORMATC/'  ENTER  HEAT  TRANSFER  COEFFICIENT  OF  ', 

*  •GUTER  LIQUID  SURFACE:') 
READ(IN,2)  HS 

WRITEdOTf  10) 

10  FORMATC/'     ENTER    THERMAL    CONDUCTIVITY    OF    FINS     AND     ', 
*•     TUBE:1) 

READ(IN,2)    CTHM 

WRITE(IOTfll) 

11  FORMATC/'     ENTER    THERMAL    CONDUCTIVITY    OF     FLUID    ■     , 

*  ■     SURROUNDING    TJBE    AND    FINS:') 
READ(IN,2)    CTHL 

WRITE<I0T,12) 

12  FORMATC/'     ENTER    TEMPERATURE    OF    FLUID    INSIDE    TUBE:*) 
READ(IN,2)    TIN 

WRITECIQT.13) 

13  FORMATC/'     ENTER    TEMPERATURE    EXTERNAL    TO    ', 
*'     SURROUNDING    FLUID:') 

RE AD (IN, 2  J    TS 
WRITECI0T,14) 

14  F0RMAT(/'     ENTER    THE    NUMBER    OF    COLUMNS    DESIRED    ', 
*•     FOR    YOUR    MODEL:      ') 

READ(IN»15)    N 

15  FORMAT!  BN,  1 2) 

WRITECIOT,  16) 

16  FORMAT!/'     ENTER    THE   NUMBER    OF    ROWS    DESIRED    FOR    ', 
*•     THE    TUBE:' ) 

READ( IN,15)U 
WRITEdOT, 99) 
99         FORMATC/1     ENTER    THE    NUMBER    OF    COLUMNS    DESIRED    FOR    ' 
*•     THE    FIN:     •) 
READ(IN,15)ZW 
C 

CO=INT(300.0/FLOAT(N) ) 
C  INITIALIZE    THE    P    AMD    CK    ARRAYS 

C 


DO   25    I=1,C0*N 

DO   25   J=lt4 

P(  I, J)  =  0 

CK(I  ,  J)=0.0 

25 

CONTINUE 

C 

c 

VALUES    COMMON    TO 

c 

L21=C0*N 

L22=2 

L23=0 
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L24=0 

L25=0 

L26=0 

L2  7=0 

L28=0 
C 

L31=300 

L32=50 

L33=6 

L34=2 

L3  5=4 

L36=6 
C 

FL41=.01 

FL42=. 66667 

IL43  =  12 

FL44  =  .8 

FL45=100. 
C 

LI  NT  1=0 

LINT2=0 

L  I  NT  3=0 
C 

C  COMPUTE    INTERNODE    CONNECTIONS 

C 

DO    55    I  =  l,CO*N 

X  =  N*(C0-1)  +  1 

IF(AMOD(FLOAT(  1-1 )  ,FLOAr(N)  J.EQ.OJ     GO    TO    30 

IF(AMOD<FLGAT(I)tFLOAT(NI).EQ.Q)    G3    TO   40 

IF(I.LE.N)    G3    TO    42 

IF(I.GE.X)    GO    TO    48 
C 
C 


P(Itl)=(I-NJ*10+l 
P(  It 2)  =  CI-1)*10+1 

P(I,3)=(I+1)*10+1 

P(  I,4)=U  +  N)*10+1 

GO   TO    55 
C 
C 
30  IFU.LE.N)    G3    TO    45 

IF(I.GE.X)    GO    TO    49 

P(I,2)=(H-1)*10+1 
35  P(Itl)=(I-N)*10+l 

P(  I,3)=(I+NJ*10+1 

GO    TO    55 
C 
C 
40  IF(I.GE.X)    G3    TO    51 

IF(I.LE.N)    G3    TO    47 

P<  I,2)=(I-1)*10+1 

GO   TO    35 
C 
C 
42  P(I,1)=(I-1)*10+1 

P(  It2)=(I  +  l)*10+l 

P<It3)=(I+N)*10+l 

P(I,4)=3014 

GO    TO    55 
C 
C 

45  P(  I, 1)=(I*20)  +  1 

46  P(If2)=lI+N)*10+l 
P(  I,3)=3014 

GO    TO    55 
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c 
c 

47  PCIt 1J=( I-l)*10+l 
GO    TO    46 

c 

48  P(  If  1)=(I-N)*10+1 
P(Ii 2)=(I-1)*10+1 

P{  If3)=(I+l)*10+l 

P( It4)=3024 

GO    TO    55 
C 
C 

49  P<  I,2)=(I-1)*10+1 

50  P(Itl)  =  U-NJ*10«-l 
P(  I,3)  =  3024 

GO    TO    55 
C 
C 

51  P(  I,2)=(I-1)*1(H1 
GO    TO    50 

C 

c 

55  CONTINUE 

C 

C  CALCULATE    HEIGHT     AND    WIDTH    INCREMENTS 

C 

C  ZW=INT[FL0AT(N)/3.0) 

ZW1=N-ZW-1 
C  ZW2=INTUFL0AT(C0)-10.0)  /2.0) 

ZW2=2 

ZW3=C0-(U+1)-ZW2 

WINC1=D/FL0AT( ZWJ 

WINC2=C/FL0AT( ZW1  ) 

HINC1=(R0-RI)/FL0AT(U) 

HINC2=B/FL0AT( ZW3 ) 

HINC3=A/FL0AT( ZW2 ) 

EL=( WINC1+WINC2J/2.0 
C 

C  COMPUTE    THE    RADIUS    OF    EACH    ROW    OF    NODES 

C 

DO    1000    1=1, DO 

IFd.GT.U    GO    TO    100 

RADH=RI*0.5*HINC1 

RADL=RI 

GO    TO    150 
C 
100         IF (I .GE.(U+1) )    GO    TO    105 

RADH=RI  +  (FLOATU)-0.5)*HINC1 

RADL=RI  +  < FLOAT ( I ) -1 .5) *H I NC  1 

GQ   TO    150 
C 
105         IF(I.NE.(U+1) )    GO    TO    110 

RADH=RO+HINC2/2.0 

RADL=RH-(FL0ATU)-1.5)*HINC1 

GO    TO    150 
C 
110         IF(I.GE.(11+ZW3JJ     GO    TO    115 

R  ADH=RO  +  <  FLOAT  m-(FLOAF  (UJ+0.51  )*H1NC2 

RADL=RO+iFLOAT<I I -(FLOAT (U) +1.5))*HINC2 

GO    TO    150 
C 
115         IF(I.NE.(U+1+ZW3) J    GO    T3     120 

RADH  =  R0  +  B  +  HINC3/2  .0 

RADL=R0+6-HINC2/2 .0 
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GO    TO    150 

C 

120         IF(I.EQ.CO)    GO    TO     125 

RADH=RQ*-B+(FLOAT(  I)-(FL3AT(  U)  +0.  5  )  -FLOAT  (  ZW3  )  )*HINC3 

RADL=RO+B+(FLOAT< I )-(FL3AT( U) +1 . 5) -FLO AT ( ZH3 ) )*HI NC3 

GO    TO    150 
C 
12  5         RAOH=RO+B+A 

RADL=R0+b+A-HINC3/2.0 
C 
C 
C 

150         ALINE=RADH**2-RADL**2 
C 

C  COMPUTE     INTER-NODE    CONDUCTANCE    \/ALJES 

C 
C 
18  5         IF(I.GT.l)    GO    TO     187 

HTC=HI 

RD=RI 

GO    TO    189 
187         iF  (I.NE.CO)    GO    TO     190 

HTC=HS 

RD=RO+A+E 

189  Al=HTC*WINCl*RD/2 .0 
A2=A 1*2.0 
A3=A2*EL/KINC1 
A4=HTC*WINC2*RD 
A5=A4/2.0 

C 
C 

190  IF(I.GE.(U+2) )    GO    TO    203 
C 

VAL4=RADL*WINC1*CTHM/HIMC1 

GO    TO    300 
200         IFU.GE.<(U+2)+ZW3))    GO    TO    250 

VAL4  =  RADL*WINC1*CTHM/HI\IC2 

GO    TO   300 
250         VAL4=CTHL*RADL*WINC1/HINC3 
300         VAL2=VAL4/2.0 
C 
C 
C 

c 

I F ( I .GE.(U+1J )    GO     TO    310 

VAL14=CTHM*RADH*W INC1/HINC1 

GO    TO   400 
310         IF(I.LE.(U+ZW3))     GO    TO    350 

VAL14=CTHL*RADH*H INC1/HINC3 

GO    TO   400 
350         VAL14=CTHM*RADH*WINC1/HINC2 
400         VAL12=VAL14/2.0 
C 
C 
C 

IFU.GE.<U+2)  )    GO    TO    413 

VAL6=CTHM*RADL*EL /HINC1 

GO   TO    5  0C 
410         IFU.GE.UU+2)+ZW3)>    GO    TO    450 

VAL6=((CTHM*WINC1+CTHL*W  INC2)*RADL)/  (2.0*HIMC2) 

GO    TO    500 
450         VAL6=CTHL*RADL*EL/HINC3 
C 
C 
C 
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500         IF (I .GE. ( U+l ) )    GO     TO    613 

VAL16=CTHM*EL*RADH/HINC1 

GO    TO   650 
610         IF(I.LE.(U+ZW3)J    GO    TO    620 

VAL16=CTHL*EL*RADH/HINC3 

GO    TO    650 
62  0         VAL16  =  (  <CTHM*WINC  I  +CTHL*  M  INC2  l*\  ADH  )/  (  2.0*HI  NC2  J 
C 
C 
C 
650         IF(I.GE.(U+2H    GO    TO    670 

VAL8=RADL*WINC2*CTHM/HINCl 

GO    TO    700 
670  IF ( I .GE.((U+2)+ZW3  ))    GO    TO    680 

VAL8=CTHL*RADL*WINC2/HINC2 

GO    TO    700 
680         VAL8=CTHL*RADL*WINC2/HINC3 
700         VAL10=VAL8/2.0 
C 
C 
C 

IF(I  .GE.(U-HJ)    GO    TO    730 

VAL18=CTHM*RADH*WINC2/HINC1 

GO    TO    750 
730         IF(I.LE.(U+ZW3))     GO    TO    740 

VAL18=CTHL*RADH*WINC2/HINC3 

GO    TO    750 
740         VAL18=CTHL*RADH*WINC2/HINC2 
750         VAL20=VAL18/2.0 
C 
C 

c 

IF(I.NE.UU+1)+ZW3)J    GO    TO    760 

VAL5=(CTHL*(RADH**2-(R0+B)**2)+CTHM*( (RO  +  B) 
***2-RADL**2))/(2.0*WINCl ) 

GO    TO    800 
760         IFU.GT.UU+U+ZW3))    GO    TO    770 

VAL5=CTHM*ALINE/( 2.0*WINC1) 

GO    TO    800 
770         VAL5=CTHL*ALlNE/( 2.0*WINC1) 
C 
C 
C 
800         IFU.GE.(U  +  1J)    GO     TO    810 

VAL7=CTHM*ALINE/<  2.0*WI^C2) 

GO    TO    845 
810         IF(I.NE.<U+1) )    GO    TO    823 

VAL7=(CTHL*(RADH**2-R0**2)+CTHM*(R0**2-RADL**2)  ) 
*/ <2.0*WINC2) 

GO    TO    845 
320         VAL7=CTHL*ALlNE/( 2.0*WINC2) 
C 
C  COUPLE    CONDUCTANCE    VALUES    TO   CONNECTION    TERMS 

C 

845         IF(I.EQ.l)    GD    TO    905 

IFU.EQ.CO)    GO    TD     912 

0=N*(I-1) 
850         DO    900    J=Q+1.Q+N 
C 
C 
860         IF(J  .EU.  (Q  +  l)  J    GO     TO    875 

IFIJ.LT«(Q+Zri+U)     GO    TO    880 

IF(J  .EQ. lg  +  Za  +  1) J      GO    TO    885 

IFCJ.NE. (Q+Ni  )     GO     TO    890 
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CK(J,1J=VAL13 

CK(J,2)=VAL7 

CK(J,3)=VAL20 

GO  TO  900 
C 
C 
880    CK(J,1)=VAL4 

CMJ,2)=VAL5 

CK(J,3)=VAL5 

CK(J,4)=VAL14 

GO  TO  900 
C 
C 
885    CK(J,1)=VAL6 

CK(J,2)=VAL5 

CK(J,3)=VAL7 

CK(J,4)=VAL16 

GO  TO  900 
C 

c 

890    CK(J,1)=VAL8 

CK(J,2)=VAL7 

CK( J,3J=VAL7 

CK(Jf4)=VAL18 

GO  TO  900 
C 
C 
875    CK(J,U=VAL2 

CK(J, 2)=VAL5 

CK(J,3)=VAL12 
C 

900  CONTINUE 

901  GO    TO    1000 
C 

c 

905  DO    910    J=1,N 
IF(J.GT.l)    GO    TO    9  06 
CK(J,  1J=VAL5 
CK(J,2)=VAL12 
CK(J,3)=A1 

GO    TO    910 
C 

c 

906  IF(J.GE.ZW+1)     GO    TO   907 
CK(J,1)=VAL5 
CK(J,2)=VAL5 

CK( J,3)=VAL14 

CMJ,4)=A2 

GO    TO    910 
C 
C 

907  IFCJ.NE.CZW+1) J    GO    TO    908 
CK(J,1)=VAL5 
CK(J,2J=VAL7 
CK(J,3)=VAL16 
CK(J,4)=A3 

GO    TO    910 
C 
C 

908  IFU.EQ.N)     GO    TO    909 
CK(J,1)=VAL7 
CK(J,2)=VAL7 
CK<J,3J=VAL18 

CK(  J,4)=A4 
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GO    TO    910 
C 
C 

909  CK(J,1)=VAL7 
CK(J,2)=VAL20 
CK(J,3)=A5 

C 

910  CONTINUE 
C 

c 

91  1         GO    TO    1000 

C 

C 

912  Q=(C0-1)*N 

DO    917    J=Q+1,Q+N 

IF(J.GT.(Q+1) )    GO    TO    913 

CK(J,1)=VAL2 

CK(J,2)=VAL5 

CK(J,3)=A1 

GO    TO   917 
C 
C 

913  IF(J .GE.(ZW+Q+l) J     GO    TO    914 
CK(J,1)=VAL4 
CK(J,2)=VAL5 
CK(J,3)=VAL5 

CK(J,4)=A2 
GO    TO    917 
C 

c 

914  IF(J.NE.(ZW+Q  +  1H     GO    TO    915 

CK(J,1J=VAL6 
CK(J,2)=VAL5 
CMJ,3)=VAL7 
CK(J,4)=A3 
GO    TO    917 
C 

c 

915  IFCJ.EQ.CQ+NJ )    GO    TO    916 
CK(J,1)=VAL8 

CK(J, 2)=VAL7 
CK(J,3)=VAL7 
CK(J,4)=A4 
GO    TO    917 
C 

c 

916  CK(J,1)=VAL10 
CK(J,2)=VAL7 
CK( Jt3i=A5 

917  CONTINUE 
1000  CONTINUE 
C 

c 

C         END    OF    CONDUCTANCE    CALCULATIONS 
WRITE(I0T,922) 

922  FORMAT!/, ■     ENTER    NAME    OF     INPUT    FILE1:) 
READ(IN,923 )FNAME 

923  F0RMAT(A12) 
C 

C 

C         CREATE    DATA    FILE    FOR    TVSSI 


C 


OPEN  (I WR,  F  I  LE=F  NAME,  STATUS^  NEW  ,FO*M=  •  F  ORMATTE  D' 
*IOSTAT=ICK) 
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WRITE(I0T,925) 

925  FORMAT(/«     WRITING    FILE    TVSSI1) 
C 

WRITE(I0T,926)C0*N,N,C0, ZW2 

926  FORMAT(/»     THIS     '.ISt'     NODE    MODEL    HAS    •,12,'     RQriS     • 
*■     AND    tfl2tt    COLUMNS.       THERE    ARE    • ,12,'    ROwS ■ ) 

WRITE(IQT,927)ZW3 

927  FORMATd     IN    THE    CONDESATE    LAYER    AND     "flZt1    ROWS    ■, 
*•     IN    THE    FIN.     THE     FIN    CONTAINS1) 

WRITE(IOT,928) ZW 

928  F0RMAT(I2,»    COLUMNS.     ■  J 
C 

C         WRITE    STATEMENTS 
C 

WRITE(IWR,943)     TITLE 

WRITE (I WR, 944  J L21 , L22,L23,L 24, L2 5, L26, L27, L28 

WRITE  CI WRt950 ) LINT1,LINT2 ,L INT3 

WRITE(IWR,946)L31,L32,L3  3,L34,L35,L36 

WRITE(IWR,947)FL4  1,FL42, I L43,FL44, FL45 

WRITE(IWR,943)TI,  TS 
C 

DO    1100    I=1,C0*N 

IF    (Pdf4i.EQ.0J    GO    TO    1011 

IF(CK(I,1).GE.  100  .O.OR.SKd  ,2  )  .GE.  100.  0)  GO    TO    1005 

IF(CK(I,3).GE.100 .O.OR.CKd ,4) .GE. 100 . 0)G0    TO    1005 

WRITE (IWR, 948 JP( I ,1),P< I ,2) ,P( 1,3) ,P( 1,4) 

WRITE (IWR,960) CM  I , 1) ,CK ( I , 2) ,CK (I ,3),CK(I,4) 

GO    TO    1100 
10  05      WRITE(IWR,948)P( I  ,  l),P(I  ,2)  ,P(I, 3)  ,P( 1,4) 

WRITE(IWR,942)CK( I , 1) ,CK ( I , 2) ,CK( I ,3 ) ,CK(I,4) 

GO    TO    1100 
1011       IF(CK(I  ,1).GE  .100.  O.OR.CKd  ,2).  GE.  100.0)     GO    TO    1015 

IF(CKd,3).GE.  103.0)    GO    TO    1015 

WRITE (IWR, 945 )P( I ,1),P(I  ,2)  ,P(I,3) 

WRITE(IWR,949)CK( I , 1 ) ,CK ( I , 2) ,CK ( I , 3 ) 

GO    TO   1100 
10  15       WRITE (IWR, 945 )P( I ,1),P(I ,2) ,P(I,3) 

WRITE (IWR, 980 )CK( 1,1),  CM I,2),CK(I,3) 
C 

1100      CONTINUE 
C 

C         FORMAT    STATEMEMTS 
C 

940  F0RMAT(1X,A79) 

941  F0RMAT(6F8.5) 

942  F0RMAT(4F8.4) 

943  FORMAT(2F8.0) 
C 

944  F0PMAT19I4) 

945  F0RMAT(9I3) 

946  F0RMAT(9I4) 

947  F0RMAT(F8.3,F8.5, I8.2F8.  1) 

948  F0RMAT(9I8) 

949  F0RMAT(F8.6,F8.5,F8.6) 

950  F0RMAT(9I4) 

960         F0RMAT(F8.6,2F8.5,F8.6) 

980         F0RMAT(3F8.4) 

C 

C 

CL0SE(IWR,I0STAT=ICK) 

IFdCK    .NE.    0)    WRITE(I0T,921)EEX 
921         FORMATt 'TROUBLE    CLOSING    PRINTER    OUTPUT    FILE*) 
C 
C 
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STOP 

END 
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A££EUCLU  fi 
Ql    COMPUTER    GENERATED    MODEL    VERIFIER    PROGRAM 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C  PROGRAM  Ql  C 

C  PROGRAM  Ql  IS  DESIGNED  TO  ACCEPT  VARIABLE  C 

C  DIMENSIONS,  INPUT  A  DATA  FILE  OF  NODE  C 

C  TEMPERATURES  GENERATED  BY  TVSSI  AND  CREATE  C 

C  AN  OUPUr  FILE  OF  HEAT  FLOW  VALUES  ACROSS  C 

C  1)  VAPOR  CONDENSATE  INTERFACE,  2)  FIN  AND  C 

C  TUBE  OUTER  PERIMETER,  AND  TUBE  INNE*  WALL  C 

C  C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c 

DIMENSION    CK(300,4),T!300) 

INTEGER    P!300,4),C0,N,X,ZW,ZW1,ZW2,ZW3,V1,V2  ,W,U, 
DOUBLE    PRECISION    CTHM ,CT HL, HI ,HS ,RO, RI , A ,B, C , D, 
*VAL4,VAL14,VAL6,VAL16,VAL8, VAL18 , VAL5 , VAL7, V ALIO , 
*A2,A3,A4,A5,WINC1 , WINC2, HINC1 ,HI NC2, HI NC3 , AL INE , 
*TIN,TS,T,X1.X2,VAL2,VAL12,VAL20,A1,RADH, RADL , EL 
CHARACTER    *79    TITLE 
CHARACTER    *79    FNAME 
C 

DATA     IOT.IN, IPR,I WR/6,5,4,8/ 
C 

OPEN{IPR,FILE=,PRN,,STATUS=,NEWl  , FQ*M  = ■ FORMATTED*  , 
*IOSTAT=ICK) 
IFUCK    .NE.    0)    WRITE(  IOT,20) 
C 

20         FORMAT*1 TROUBLE   OPENING    PRINTER    OUTPUT    FILE1) 
C 

WRITEUOT,  17) 

17  FORMAT!/'     ENTER    THE    TITLE    OF    THIS    STUDY,    79     ', 

*  '    COLUMNS    MAX.       THIS    TITLE    WILL    APPEAR1, /,•     AT', 

*  ■    THE    TUP    OF    EVERY   PRINTED    PAGE    OF    OUTPUT: •) 
READ(IN,18)    TITLE 

18  F0RMAT<A79) 


C 


WRITEUOT, 1) 

FORMAT!/'     ENTER    INNER    RADIUS    OF    TUBE:*) 

READ!IN,2)    RI 

F0RMAT(BN,F10.5) 

WRITEUOT, 3) 

FORMAT!/'     ENTER    OUTER    RADIUS    OF    TUBE:1 J 

READ!IN,2)    RO 

WRITEUOT, 4) 

FORMAT!/1     ENTER    DISTANCE    FROM    OUTER    EDGE    OF 
*     'FIN    TO    LIQUID    SURFACE:'] 
READUN,2)    A 

WRITEUOT, 5) 

FORMAT!/'     ENTER    HEIGHT    OF    FIN') 

READ (IN, 2)     B 

WRITEUOT, 6) 

FORMAT!/'     ENTER    DISTANCE    BETWEEN    FINS') 

READ(IN,2)    C 
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c 


c 


c 


WRITEdOT, 7) 

7  FORMATC/'     ENTER    WIDTH    OF    FIN') 
READ(IN,2)    0 

WRITEdOT, 8) 

8  FORMATi/*     ENTER   HEAT    TRANSFER    COEFFICIENT    2F     ■, 

*  '     TUBE    INNER    SURFACE:  ■  ) 
READ<IN,2)    HI 

WRITEdOT, 9) 

9  FORMAT!/'     ENTER    HEAT    TRANSFER    COEFFICIENT    3F     ', 

*  'OUTER    LIQUID    SURFACE:1 J 
READUN,2)    HS 

WRITEdOT, 10) 

10  FORMAT*/'     ENTER    THERMAL    CONDUCTIVITY    OF    FINS    AND    ', 
*•     TUBE:*) 

READ(IN,2)    CTHM 

WRITE(I0T,11) 

11  FORMAT*/'     ENTER    THERMAL    CONDUCTIVITY    OF    FLUID    ■     , 

*  •    SURROUNDING    TUBE    AND    FINS:1) 
READ(IN,2)    CTHL 

WRITEdOT,  12) 

12  FORMAT*/'     ENTER    TEMPERATURE    OF    FLUID     INSIDE    TUBE:') 
READ(IN,2)    TIN 

WRITEdOT, 13) 

13  FORMAT(/'     ENTER    TEMPERATURE    EXTERNAL    TO    '  , 
*•     SURROUNDING    FLUID:1) 

READUN.2)    TS 
WRITEdOT, 14) 

14  FORMAT*/'     ENTER    THE    NUMBER    OF    COLUMNS    DESIRED    ', 
*•     FOR    YOUR    MODEL:     •) 

READ*IN,15)    N 

15  FORMAT(BN,  12) 

WRITEdOT,  16) 

16  FGRMAT*/'     ENTER    THE    NUMBER    OF    ROWS    DESIRED    FOR    ', 
*•     THE    TUBE:' ) 

READ( IN,15)U 

WRITEdOT, 99) 
99         FORMAT*/1     ENTER    THE    NUMBER    OF    CDLUMNS    DESIRED    FOR    • 
*•     THE    FIN:     •) 

READ( INf 15JZW 
C 

C0=INT(300.0/FL0AT(N) ) 
C         INITIALIZE    THE    P    AND    CK    ARRAYS 
C 

DO   25    I=1,CG*N 

DO    25    J=l,4 

P(  I, J)  =  0 

CK(I  ,  J)=0.0 
25  CONTINUE 

C 

C         VALUES    COMMON    TO    ALL    RUNS 
C 

L21=C0*N 

L22=2 

L23=0 
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L24=0 

L25=0 

L26=0 

L27=0 

L28=0 

c 

L31=300 

L3  2=50 

L33=6 

L34=2 

L35=4 

L36=6 

c 

FL41=.01 

FL42  =  . 66667 

I L 43 =12 

FL44=.8 

FL45=100. 

c 

LI  NT  1=0 

LI  NT  2=0 

LI  NT  3=0 

c 

c 

COMPUTE    INTERNODE 

c 

CONNECTIONS 

DO    55    I=1,CC*N 

X=N*(C0-1)+1 

IF(AMOD(FLOAT(  1-1 ), FLOAT (N)  ).EQ.O)     30    TO    30 

IF(AMOD(FLOAT( I),FLOAT(Ni J.EQ.O)    GO    TO    40 

IF(I.LE.N)    GO    TO    42 

IF(I.GE.X)    GO    TO    48 


P(I, 1)=(I-N)*10+1 
PUi2J={  1-11*10+1 
P(  It 3)=(I  +  1)*10  +  1 

P(  I,4i=(I+N)*10+l 

GO    TO    55 

C 

c 

30  IFU.LE.N)    GO    TO    45 

IFU.GE.X)    GO    TO    49 
P(I,2)=(I+1)*10+1 

35  P( I, l)=(I-Ni*10+l 

P(  I,3)=(I+N)*10+1 
GO    TO   55 

C 

C 

40  IF(I.GE.X)    GO    TO    51 

IF(I.LE.N)    GO    TO    47 
P(  If2)=(I-l)*10+l 
GO   TO   35 

C 

C 

42  P(I,1)=(I-1)*10+1 

P( It 2J=(I+1)*10+1 
P(I,3)=(I+N#*10+1 
P(If4)=3014 
GO    TO   55 

C 

C 

45  P(  I, 1)=(I*20)+1 

46  P(  I,2)=(I  +  N)*10+1 
P(  I, 3)  =  3014 

GO    TO    55 
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c 
c 

47  P(  I,  U=(  1-1  1*10+1 
GO    TO    46 

C 
C 

48  P(  I,l»=( I-N)*10+l 
P(  It  2  )=C 1-11*10*1 

P(  It 3)  =  ( 1+1 1*10+1 
P(  If4)=3024 

GO    TO    55 
C 
C 

49  P(I,2)=(I-1)*10+1 

50  P(  I,l)=(  I-N)*10+l 
P(  I,3)  =  3024 

GO   TO    55 
C 
C 

51  P( I, 2)=(I-1)*10+1 
GO    TO    50 

C 

C 

55  CONTINUE 

C 

C  CALCULATE    HEIGHT    AND    WIDTH     INCREMENTS 

C 

C  ZW=INT(FL0AT(N)/3.0) 

ZW1=N-ZW-1 
C  ZW2=INT((FL0AT(C0)-10.0)/2.0) 

ZW2=2 

ZW3=C0-(U+1)-ZW2 

WINC1=D/FL0AT(  ZWI 

WINC2=C/FL0AT( ZWI ) 

HINC1=(R0-RI)/FL0AT(U) 

HINC2=B/FL0AT(ZW3) 

HINC3=A/FL0AT( ZW2 ) 

EL=( WINCl+WlNC2i/2.0 
C 

C  COMPUTE    THE    RADIUS    OF    EACH    ROW    OF    MODES 

C 

DO    1000    1=1, CO 

IF(I.GT.l)    GO    TO     100 

RA0H=RI+0.5*HINC1 

RADL=RI 

GO    TO    150 
C 
100         IF(I.GE.(U+1) )    GO    TO    105 

RADH=RI+(FL0AT(I)-0.5)*HINC  1 

RADL=RI+(FL0AT(I)-1.5)*HINC1 

GO    TO    150 
C 
105         IF(I.NE.(U  +  1)  )    GO     TO    110 

RADH=RO+HINC2/2.0 

RADL=RI  +  <FL0ATU)-1.5)*HINC1 

GO  TO  150 
C 
110    IFU.GE.UH-ZW3))  GO  TO  115 

RADH=RO+( FLOAT (I) -(FLOAT <U) +0.5) )*HINC2 

RADL=RO+( FLOAT ( I ) - (FLOAT ( U) +1 .5 ) }*HINC2 

GO    TO    150 
C 
115         IF(I  .NE.(U+1+ZW3J  )    GO    TO     120 

RADH=R0+B+HINC3/2.0 

RADL=R0+B-HINC2/2 .0 
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GO    TO    150 
C 
120         IF(I.EQ.CO)    GO    TO     125 

RADH=RO+B+(FLOAT(I )-{ FLOAT ( U)+0. 5 ) -FLOAT ( ZW3 ))* HI NC3 

RADL  =  RO+B  +  (  FLOAT! I)-(FL3AT( U)  +  l.  5) -F LOAT ( ZK3 ) )*HINC3 

GO    TO    150 
C 
125         RADH=RO+B+A 

RADL=R0+B+A-HINC3/2.0 
C 
C 
C 

15  0         ALINE=RADH**2-RADL**2 
C 

C  COMPUTE    INTER-NODE    CONDUCTANCE    VALUES 

C 
C 
18  5         IF(I.GT.l)    GO    TO    187 

HTC=HI 

RD=RI 

GO    TO    189 
187         IFU.NE.CO)    30    TO     190 

HTC=HS 

RD=RO-t-A+B 

189  Al=HTC*WINCl*RD/2.0 
A2=A1*2.C 
A3=A2*EL/HINC1 
A4=HTC*WINC2*RD 
A5=A4/2.0 

C 
C 

190  IFU.GE.CU+2)  J    GO    TO    203 
C 

VAL4=RADL*W INC 1*C THM/HINC 1 

GO    TO    300 
200         IF(I.GE.((U+2)+ZW3J)    GO    TO    250 

VAL4=RADL*W INC 1*CTHM/ HINC2 

GO    TO   300 
250         VAL4=CTHL*RADL*WINC1/HINC3 
300         VAL2=VAL4/2.0 
C 
C 
C 

c 

IF(I.GE.(U+1) J    GO     TO    310 

VAL14=CTHM*RADH*K INC1/HINC1 

GO    TO    400 
310         IF(I.LE.(U+ZW3 ))     GO   TO    350 

VAL14=CTHL*RADH*WINC1/HINC3 

GO    TO    400 
350         VAL14=CTHM*RADH*WINC1/HINC2 
400         VAL12=VALl4/2.0 
C 
C 
C 

IF(I.GE.(U+2) )    GO    TO    413 

VAL6=CTHM*RADL*EL/HINC1 

GO    TO    500 
410         IF(I.GE.<(U+2)+ZW3J)    GO    TO    450 

VAL6  =  UCTHM*WINC1+CTHL*WINC2)*RADL)/  (2.0*HI^C2) 

GO    TO    500 
450         VAL6=CTHL*RADL*EL/HINC3 
C 
C 
C 
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500         IF(I.GE.(U  +  1)  )    GO    TO    610 

VAL16=CTHM*EL*RADH/HINC1 

GO    TO   650 
610         IFU.LE.IU+ZW3  ))     GO    TO    620 

VAL16=CTHL*EL*RADH/HINC3 

GO    TO    650 
620         VAL16  =  (  (CTHM*WINC1+CTHL*WINC2)*RADH  )/  (2.0*HINC2) 
C 
C 

c 

650         IFCI.GE.1U+2) J    GO     TO    670 

VAL8=RADL*WINC2*CTHM/HINC1 

GO    TO    700 
670         IF(I.GE.((U+2)+ZW3)J    GO    TO    680 

VAL8=CTHL*RADL*WINC2/HIMC2 

GO    TO    700 
680         VAL8=CTHL*RA0L*WINC2/HINC3 
700         VAL10=VAL8/2.0 
C 
C 
C 

IFCI.GE.(U+1) )    GO    TO    730 

V AL1 8=C  THM*RADH* W INC2/HI  NCI 

GO    TO    750 
730         IF  (I.LE.CU  +  ZW3 J)     GO    TO    740 

VAL18=CTHL*RADH*W INC2/HI  NC3 

GO    TO    750 
740         VAL18=CTHL*RADH*WINC2/HINC2 
750         VAL20=VAL18/2.0 
C 
C 
C 

IF(I.NE.UU+U+ZW3))    GO    TO    760 

VAL5=(CTHL*(RADH**2-(R0+B)**2)+CTHM*( (RO+B) 
***2-RAUL**2))/ (2.0*WINC1  ) 

GO    TO    800 
760         IF(I.GT.((U+1)+ZU3))    GO    TO    770 

VAL5=CTHM*ALINE/(2.0*WINC1) 

GO    TO    800 
770         VAL5=CTHL*ALINE/(2.0*WINC1J 
C 
C 
C 
800         IF(I.GE.(U  +  m    GO     TO    810 

VAL7  =  CTHM*ALINE/(  2.0*WI^C2) 

GO    TO    845 
810         IFU.NE.lU  +  1)  )    GO    TO    82Q 

VAL7=<CTHL*(RADH**2-R0**2  )+CTHM*( R0**2-R ADL**2)  ) 
*/ (2.0*WINC2) 

GO    TO    845 
820         VAL7=CTHL*ALINE/(2.0*WINC2) 
C 

C  COUPLE    CONDUCTANCE    VALUES    TO   CONNECTION    TERMS 

C 
845         IFCI.EQ.1J     GO    TO    905 

IF(I.EQ.CO)    GO    TO    912 

Q  =  N*U-1) 
850         DO    900    J=Q+1,Q+N 
C 

c 

860         IF(J.EQ.(Q+1)  )    GO  TO    375 

IF(J.LT.(Q+ZH+1)J  GO    TO    880 

IF(J.EQ.(Q+ZM+l))  GO    TO    885 

IF (J .NE. (Q+NJ )     GO  TO    893 
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c 
c 

880 


C 
C 
885 


C 
C 
890 


C 
C 
875 


C 

900 
901 
C 

c 

905 


C 
C 

906 


C 
C 
907 


C 

c 

908 


CK(J,1)=VAL1D 
CKU,2)=VAL7 
CMJ,3)=VAL20 
GO    TO    900 


CK( J,1)=VAL4 
CK(J,2)=VAL5 
CK(J,3)=VAL5 
CK(J,4)=VAL14 
GO    TO    900 


CK(J,1)=VAL6 
CK(J,2)=VAL5 
CK(J,3)=VAL7 
CK(J,4J=VAL15 

GO    TO    900 


CK(J,1J=VAL8 
CMJ,2)=VAL7 
CK(J,3)=VAL7 
CM  J,4)=VAL18 

GO    TO    900 


CK(J,1)=VAL2 
CK(J, 2)=VAL5 
CK(J,3)=VAL12 

CONTINUE 
GO    TO    1000 


DO   910    J=1,N 

IF(J.GT.l)    GO    TO    906 

CK(J,1)=VAL5 

CK<J,2)=VAL12 

CK(J,3)=A1 

GO    TO    910 


IFU.GE.ZW+1)     GO    TO   907 

CK(J,1)=VAL5 

CK( J,2J=VAL5 

CK<J,3)=VAL14 

CK( J,4)=A2 

GO    TO    910 


IF(J.NE.(ZW  +  1)  ) 

CKU,  1)=VAL5 

CK(J,2)=VAL7 

CK(J,3)=VAL16 

CK{Jt4)=A3 

GO    TO    910 


GO   TO    908 


IFU.EQ.N)    GO 

CK(J,1)=VAL7 

CK(J,2)=VAL7 

CK(J,3)=VAL18 

CK(J,4)=A4 


TO    909 
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GO  TO  910 
C 
C 

909  CK(J,1)=VAL7 
CK(J,2)=VAL20 
CK(J,3)=A5 

C 

910  CONTINUE 
C 

C 

91  1    GO  TO  1000 

C 

c 

912  Q=(C0-1)*N 

DO    917    J=Q+1.Q+N 
IFU.GT.(tf+lJ  )     GO    TO    913 
CK(J,1J=VAL2 
CK(J,2)=VAL5 
CK(J,3)=A1 
GO    TO    917 
C 

c 

913  IF(J.GE.IZW+Q+1))     GO    TO    914 
CK(J,1)=VAL4 

CK( J,2)=VAL5 

CK(J,3)=VAL5 

CK(J,4)=A2 

GO    TO    917 
C 
C 

914  IF(J.NE.(ZW+3  +  UJ     GO    TO    915 
CK(J,1)=VAL6 
CKU,2)=VAL5 
CK(J,3)=VAL7 

CK(J,4)=A3 

GO    TO    917 
C 
C 

915  IF(J.EQ.(Q+N) )    GO    TO    916 
CK(J,U=VAL8 
CK(J,2)=VAL7 
CK(J,3)=VAL7 

CK( J,4)=A4 

GO    TO    917 
C 
C 

916  CKU,1)  =  VAL10 
CK<J,2)=VAL7 
CM  J,3)=A5 

917  CONTINUE 
1000  CONTINUE 
C 

C 

C         ENO    OF    CONDUCTANCE    CALCULATIONS 
WRITE(I0T,922) 

922  FORMAT!/,1     ENTER    NAME    OF     INPUT    FILE':) 
READ(IN,923)FNAME 

923  F0RMAT(A12) 
C 

c 

C    CREATE  DATA  FILE  FOR  TVSSI 


C 


OPEN (IWR,FILE  =  FNAME,STATUS=•NEW,  , F OR M=« FORMATTED* 
*IOSTAT=ICK) 
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C  READ    IN    NOOE    TEMPERATURES    GENERATED    BY    TVSSI 

C 


C 

924         CONTINUE 


DO   924   J  =  l,25 

READ (16 ,1045)  (T(I  )  ,  I=( ( J-l) *12+1 ) , J* 12 ) 


WRITEUQT.925) 

925  FORMAT</«     WRITING     FILE    TVSSI') 
C 

WRITc(I0T,926)C0*N,C0,N, ZW2 

926  FORMATC/'     THIS     Sn,'     NODE    MODEL    HAS     l  ilZ*9     ROrf 
*•     AND1     ,12,'     COLUMNS.       THERE    ARE    ',12,'    ROWS') 

WRITE(IGT,927)ZW3 

927  FORMATC     IN    THE   CONDESATE    LAYER    AND     ',12, 
*•     IN    THE    FIN.       THE    FIN    CONTAINS') 

WRITE(I0T,928)ZW 

928  FORMATU2,'    COLUMNS.    ') 
C 

C  COMPUTE    HEAT    FLOW    ACROSS    DESIGNATED    SURFACES 

C 

QFL=0.0 

QFS=0.0 

QFU=0.0 

QU=0.0 

QO=0.0 
C 

c 

DO    2000    1=1, !i 
949         IF(I.GT.(N-ZW-1) )     GO    TO    955 

V1=(U+1)*N+1+ZW+I 

V2=( U*N)+1+ZW+I 

IF(V2.EQ. (U+1)*N)     GO   TO    950 

QFL=OFL+(T< VI )-T( V2) ) *CK( V2 ,4) 

GO    TO    955 
9  50         QFL=QFL+(T(V1 )-T( V2) ) *CK ( V2 ,3 ) 
955         IFd.GT.ZW3)    GO    TO    960 

W=N*(U+I )+2+ZW 

QFS=QFS  +  (T(W)-TiW-l))*CK(W,  3) 
960         IF(I.GT.(ZW+1) )    GO    TO    970 

Xl=( U+ZW3)*N+I 

X2=( U+ZW3+1)*N-H 

IF(I.EQ.l)    GO    TO    965 

QFL=QFU  +  (T(X2)-T(  XI) )*CK(X1 ,4) 

WRITE(6,9000)QFU,  T (X2 ) ,T ( XI ) , CK( XI,  4 ) , XI ,X2 

GO    TO    970 
96  5         QFU=0FU+(T(X2 )-T( XI) ) *CK ( XI , 3 ) 

WRITt(6,9000)QFU,T(X2) ,T ( XI  I , CK( XI, 3 ) , XI , X2 
9000       FORMAT(1X,4F10.5,2I3) 
970         IFd.GT.l  )G0    TO    980 
975         Q0=Q0+mi)-TIN)*CK(I,3) 

QU=QU+(TS-T  <H-(N*(C0-1)  )  )  )*CK  ( (  I  +  (  N*  (CO-  1)  )  )  ,  3) 

GO    TO    2000 
980         IF(I.EQ.N)    GO    TO    975 

Q0=Q0+(T(I)-TIN)*CK(I,4) 

QU=QU+(  TS-T(I+(N*(C0-1))  )  )*CK  (  ( I  -M  N*(CG-1)  )  )  ,4) 
2000       CONTINUE 

QT1=QFU+QFS 

QTC=QT1+QFL 
C         WRITE    STATEMENTS 

WRITE(IWR,1030) 

WRITEUWR,1031)RI , RO, B 

WRITE(IWR,1032) 

WRITE(1WR,1033)    CD, A 
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WRITE!IWR,1034) 
WRITE!IWR,1035)TIN,TS 
WRITEUWR,1036) 

WRITEUWR,103  7)CTHM,CTHL,HI  ,HS 
WRITEUWR,  1047JQFU 
WRITEUWR, 1048>QFS 
WRITEUWR, 1040JQT  1 
WRITE ( I WRf 1041) QF L 
WRITEUWR, 1042)QT0 
WRITEUWR,  1043  )QU 
WRITEUWR, 1044)Q0 
C 

1030  FORMAT!/'     PIPE    INNER    RADIUS  PIPE    OUTER    RADIUS  ', 
*•     FIN    HEIGHT          •) 

1031  F0RMAT(4X,F8.5,»     M  '  ,1 IX,  F  8.  5,  •    M«  ,  5X.  F10.5,  ■     M"  ) 

1032  FORMAT!/'     FINJ    SPACING  FIN    WIDTH  ', 
*■     CONDENSATE    HEIGHT       •) 

1033  FORMATUX,F8.5,'     M  •  ,10X,  F8.  5,  •    M',7X,F9.8,'     MM 

1034  FORMAT!/'  TEMPERATURE  INSIDE  PIPE      TEMPERATURE  ', 
*•  OUTSIDE  CONDENSATE  LAYER  •) 

1035  F0RMAT(7X,F4.1,«     C  •  ,34X  ,  F4.  1,  •    C) 

1036  FORMAT!/1     K    METAL  K    LIQUID  H    INSIDE    PIPE    ', 
*•     H    IN   CONDENSATE    LAYER1) 

1037  FORMAT!F5.0,'     W/M •  ,4X , F5 . 1 ,  ■    W/M • ,3X , F 10 .0, •     W/M**2'  , 
*    5X,  10.0, •     W/M**2') 

1040  FORMAT!/'     Q    IN    ON     FIN    UPPER    AND    SIDE    SURFACES  ', 
*F9.3, •    WATTS    ' ) 

1041  FORMAT!/'     Q    IN    ON    PIPE    SURFACE  ',15X,F9.3,'     WATTS') 

1042  FORMAT!/'     TOTAL    Q     IN    THROUGH    PIPE    AVJD    FIN    SURFACE  ', 
*F5.3,  •    WATTS    •  ) 

1043  FORMAT!/'     Q    IN    THROUGH    CONDENSATE    ',12X,F9.3,'    WATTS') 

1044  FORMAT!/'     Q    OUT    THROUGH    PIPE    INNER    SURFACE       ',3X,F9.3, 
*•     WATTS     ') 

1045  F0RMATU2F6.2  ) 

1047  FORMAT!/'     Q    IN    THROUGH    JPPER    FIN    SURF ACE • , 9X , F9 . 3 , 
*•     WATTS    •) 

1048  FORMAT!/'    Q    IN    THROUGH    FIN    SIDE    SURFACE    ',9X,F9.3, 


C 

C 


*•     WATTS    •) 


CLOSE!IWR,IOSTAT= ICK) 

IF1ICK    .NE.    0)     WRITE! IOT, 1051JEEX 
1051       FCRMAT!'TROUBLE    CLOSING    PRINTER    OUTPUT    FILE') 
C 
C 

STOP 

END 
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